Theory of resonant Raman scattering of tetrahedral amorphous carbon 



O 

o 



Mickael Profeta and Francesco Mauri 
Laboratoire de Mineralogie-Cristallographie de Paris, 
Universite Pierre et Marie Curie, 4 Place Jussieu, 75252, Paris, Cedex 05, France 

(February 1, 2008) 

We present a practical method to compute the vibrational resonant Raman spectra in solids with 
delocalized excitations. We apply this approach to the study of tetrahedrous amorphous carbon. 
' We determine the vibrational eigenmodes and eigenvalues using density functional theory in the 

local density approximation, and the Raman intensities using a tight binding approximation. The 
computed spectra are in good agreement with the experimental ones measured with visible and UV 
lasers. We analyze the Raman spectra in terms of vibrational modes of microscopic units. We show 
that, at any frequency, the spectra are dominated by the stretching vibrations. We identify a very 
rapid inversion in the relative Raman intensities of the sp 2 and sp 3 carbon sites with the frequency 
of the incident laser beam. In particular, the spectra are dominated by sp 2 atoms below 4 eV and 
by sp 3 atoms above 6 eV. 
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I. INTRODUCTION 



Tetrahedral amorphous carbon is a very studied material because of its physical properties, which are close to those 
of diamond. It is composed of carbon atoms with sp 3 and sp 2 hybridization, with a fraction of sp 3 sites larger than 
40 %. The properties of tetrahedral amorphous carbon are very dependent on its microscopic structure. Thus, the 
knowledge of the microscopic structure can be used to improve the synthesis of the material. Raman spectroscopy 
have often been applied to the study of tetrahedral amorphous carbon.H However visible Raman spectroscopy is much 
more sensitive to sp 2 than to sp^|-carbons.0 Recently, it has been shown that it is possible to probe the vibrations 
of sp 3 carbons with UV Raman,tHa but few UV spectra have been presented in literature up to now. In absence of 
further experimental data, a theoretical study can elucidate the detailed dependence of the Raman spectra on the 
laser frequency. 

A theory top-compute the Raman spectra under non-resonant conditions is well established. Indeed, the Placzek 
approximation^ links the scattered Raman intensity with the electronic susceptibility of the sample. However, for 
£f~) , resonant Raman the Placzek approximation is not justified. In this work, we develop a practical method for the 
calculation of resonant Raman spectra in solids with delocalized excitations. We apply this theory to amorphous 
carbon, using density functional theory in the local density approximation to compute the vibrational modes and 
frequencies and a tight binding approximation to obtain the resonant Raman intensities. We compare our results to 
the experimental visible and UV Raman spectra and we analyze our theoretical spectra in terms of the motion of 
microscopic units. 



II. THEORETICAL RESULTS 



In Raman spectroscopy a sample of matter is irradiated by a monochromatic laser beam of pulsation cl> l , and the 
intensity of the scattered light I(lo) is measured, as a function of the pulsation u. The intensity I(ui) has a main elastic 
peak at the incident light pulsation, lu — wl, and a smaller inelastic contribution at oj ^ lul associated to energy 
transfers between the light and the sample. In this paper, we will focus on first-order vibrational Stokes spectroscopy, 
' in which the light excites a single vibrational mode. 

In order to describe the vibrational Raman scattering and to fix the notation, we need to introduce the Born- 
Oppenheimer (BO) approximation. Within the BO approximation the eigenstates of the Hamiltonian with eigenvalues 
Et can be written as \ipe,n) ® l^e.fc), where the kets IV'e.R.) and \v e ,k) are defined on the Hilbert space of the electronic 
and nuclear coordinates, respectively. The kets IV'e.R.) are the eigenvectors with eigenvalues E e (R) of the electronic 
Hamiltonian iJ ei (R), which depends parametrically on the positions R, of the N nuclei (R represents a vector of 3N 
coordinates). For every e, each E e (R) defines a BO potential energy surface. Finally the wave-functions \v e ,k) satisfy 
the eigenvalue equation: 

[T + E e (R)]\v e , k )=E%\v e , k ), (1) 
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where T is the kinetic energy of the nuclei. It is possible to define a vibrational energy as: 



et = E%-E e (R?), (2) 

where R° q is the equilibrium position of the nuclei on the BO surface E e (R) . 

In vibrational Raman spectroscopy, before the scattering, the system is on the ground state BO surface, £o(R), 
and in the vibrational state e°, after the scattering the system is still on the ground state BO surface but in a different 
vibrational state e°. In first-order Stokes Raman the harmonic approximation is assumed and the final state differs 
from the initial one by the creation of a phonon with energy lo s .0 

In non-resonant Raman scattering the energy of incident light is very far from any electronic excitation. Under 
this condition, a simple and well established theory to compute the scattering intensity I(uS), have been developed by 
Placzek.LrO Within the Placzek's approximation the intensity of scattered light is expressed in terms of the electronic 
susceptibility at the laser frequency Xu l u(wl): 



/(ul, u, lo) (X 



<9Xu lU (wl) 2 



dQ s 



((n s ) + 1)S(lo-uj l +uj s ), (3) 



where Ul and u are the polarization of the incident and scattered light, (n s ) — [exp(uj s /(kBT)) — is the average 
thermal occupation of the vibrational state s, and the derivative of the susceptibility with respect to the vibrational 
normal mode Q s is: 

dQs ~ ^ ^ dm qj ~ dH q 1 ' 

J a— x,y,z J 

Here qj' s is a component of the unit vector of normal mode Q s for vibration s on the atom J with normalization 
q s • q s = 5 S . S >. The Placzek's approximation have been used successfully for the first principles computation of the 
non-resonant vibrational Raman intensity, (see e.g. Refs. |||). 

In resonant Raman spectroscopy the energy of incident light is in resonance with an electronic transition. Under 
this condition the Placzek's approximation is, in principle, not justified. However Eq.(Q) has also been applied to 
resonant Raman scattering in solidal3 in order to interpret experimental measurements on semiconductors. 

In this paper, we will derive an expression for the scattering intensity in solids which is valid both in the non- 
resonant and resonant case. This derivation will be used in the following section to analyze the Raman spectra of 
amorphous carbon. 

We start our derivation from the generaLexpression for the vibrational Raman intensity, which can be obtained 
from theory of second quantization of light :Eil 

7(u L> u, lj) oc p{i)^\Aj?^ ~ ^ - e° + e° f ), (5) 



with 



.k 



K*l/* L (R)K fc > K*l/r(R)K/> 



S e (R c e q ) + e e k - £ (Ro q ) + e? - t^L - <7 



S e (R c e q ) + e e k - £ (R°o q ) + e° i +u + i 1 



(6) 



Here, p(i) is the probability to find the system in the initial state i, 7 is a small real number added for calculation in 
order to treat correctly the poles of the expression and 

/^(R) = 47^o,R|u-D(r,R)|^,R), (7) 
V V 



where D is the ionic and electronic dipole moment, V is the sample volume, and r is the position vectoLfor the 
electrons. Placzek's expression, Eq. (||), has been derived in the non resonant case starting from Eq. (^IxlI Notice 
that Eq. (||) is much more complicated to be evaluated than the Placzek's expression, Eq. (||), since it requires the 
knowledge of the vibrational cigenfunctions, \v e ,k), on the excited BO surfaces. nn 

Before presenting our derivation we recall, for comparison, that of the Placzek's approximation.QEl In the Placzek 
derivation Eq. (ph is approximated by: 
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^\ E e (n) - E (Il) - w h - i-y 



A,/ 



E e (R) - E (R) + cj + vy 
Using the completeness relation |f e ,fe)(i / e,ife| = 1, Eq. @ becomes: 



(8) 



, = v / K t |/u L (R)/u*(R)K/) WoMWZWKf) \ (q) 

lJ 2f \E e (K) - E (K) - cu L - i 7 E e (R) - E Q (R) + uj + ij J ' [> 

which finally leads to Eq. (|J). In order to obtain Eq. (0), we drop the dependence of the denominators of Eq. @ on 
the excited vibrational states \y e ,k)i by replacing [i? e (R° q ) + e| — i? (Rg q ) + e°] with {E e (Tl) — E (R)]- In this way, 
we neglect in the denominators an energy term 

A* = E e (fl?) - E e (R) + el- [Eo{B?) E (R) + e% (10) 

This approximation is justified if E e (R) — Eo(R) — lo^ ^> Ag. Since we are using the completeness relation, the above 
inequality should, in principle, held for every vibrational states k, and not just for the one phonon excitations. In 
particular, the relation must be verified for excitations, and energy differences A*, involving an arbitrary number of 
phonons and even for unbounded states in the energy continuum. This condition can obviously never be satisfied. 
The use of Placzek approximation relies on the hope that the inclusion of excitations involving a limited amount of 
phonons are sufficient to fulfill the completeness relation. When the laser frequency is close to an electronic transition, 
the Placzek's approximation could not be applied in general since, in this case, the dependence on Ag could not be 
ignored and the completeness relation could not be used. 

Here, to compute the Raman spectra in a periodic solid, we simplify Eq. (^) using an approximation different from 
Placzek's. In a periodic solid with dispersive bands, the electronic excitations are generally delocalized. In this case, 
the excited BO surfaces are obtained from the ground-state one by a constant vertical displacement, and we can 
assume that \v e ,k) — l^o.fc), Re q — Ro^ anc ^ e t — € k- This leads to: 



, , v ( Kii/u L (R)Kfc) K fc |/r(R)K/ 

A/(UL,U)=^^ } _ 



Ki|/°(R)K fc ) K*|/£(R)N,/) 

£ e (R° q )-£o(Ro q )+e°-e°+ W + z 7 



(11) 



We expand the expression of /„(R) at first order around Rq 01 , /„(R) — /u(Ro q ) + (R — Ro q ) ■ ^/u(Ro q )/^R an d we 
extract all the terms of degree one in (R — Rg q ) in Eq. ([n]). The term R — Rp q can be written using the operators 
of creation and annihilation of one phonon of energy u> s . If we extract Stokes terms, we can write the amplitude as a 
function of a vibration s, 



7(u L ,u,w) oc ^ p(n s )uj 4 \A s \ 2 5(uj - uj-l + ui s ), (12) 



where n s is the occupation number of the mode s, u> s = — e", and 

^ in s +T( q s -(^) R ^r(R°o q ) 



S 4" V \E e (B^) - £o(R° q ) +u s -u h -i 1 + £ e (R° q ) - E Q (K q ) - c*, - 1 7 1 dj 

q s ■ (§g) R -/-(R e q ) ^(R° q )q s • (^)r- 



S e (R^ q ) - £ (Ro q ) + w L + »7 £ e (Ro q ) - ^o(Ro q ) + + *7 ~ u s 

Finally, we neglect the dependence of the denominator over lu s . Indeed, we expect that A s has the same behavior with 
respect to as the electronic susceptibility, x(o>l). In a solid with delocalized excitations, away from the van Hove 
singularities, x(wl) is a smooth function such that xi^L ± &s) — x(^l)- Therefore, we expect that in a solid with 
delocalized excitations we can drop the term ±w s in the denominator of Eq. (|l3|). Notice that the validity requirements 
of our approximation are less stringent than those of the Placzek's approximation. Our approximation neglects the 
dependence of the denominators on the one phonon excitation energies u> S) whereas the Placzek's approximation 
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neglects the multi-phonons excitation energies A^. Therefore our approximation requires that x( w ) is on the 
energy scale of one phonon excitations whereas the Placzek's approximation requires that x(w) is flat on the energy 
scale of multiple phonon excitations. 



2uj 



1 9k(R,R',w l ) 



<9R 



(14) 



where 



(R,R',o; L )=^ 



/« (R)/r(R) 



/°(R)/£(R) 



E e (R>) - E {H') - o>l - *7 Se(R') - So(R') + w L + »7 - 



(15) 



Eqs. (|l^), (|l4|), and (|l5|) are the working expressions that we will use in the next section to compute the Raman 
intensity. 

It is interesting to compare our final expression with that of the Placzek's approximation. To this purpose we 
notice that the electronic susceptibility \ can be obtained as a special case of the function k, indeed Xu L ,u(R, wl) = 
k ULi11 (R,R,u)l). Thus, the Placzek's expression differs from the expression we derived just by an additional partial 
derivative. Indeed, if we assume that 



A 



Placzck 



2lu s 



<9k(R,R',w l ) 9k(R,R',u l ) 



<9R 



9R 



(16) 



R=R'=R° q 



we obtain Eq. (g) from Eq. |[|). 



III. CALCULATION ON AMORPHOUS CARBON 



We compute the vibrational spectra and the Raman intensities of tetrahedral amorphous carbon using a model 
generated in Ref. |l2| by tight binding (TB) molecular dynamics. The model contains 64 atoms in a periodic cubic 
cell of 7.57 A. Among the 64 atoms, 18 atoms are three-fold coordinated (sp 2 hybridized) and the others are four-fold 
coordinated (sp 3 hybridized). In this model all sp 2 atoms are inserted into a double bonded link, i.e. there are 
no electronic defects associated to dangliHS-» orbitals. This kind of defects has been identified in other computed 
generated models of amorphous carbon.li3iiJ Fig. [I] shows a ball and stick model of the sp 2 carbon atoms in the 
cell. The sp 2 atoms are arranged in four pairs and two chains of four and six atoms respectively. No aromatic ring 
are present in the sample. We compute the dynamical matrix-and the vibrational eigenvalues and eigenvectors using, 
density functional theory in the local density approximationJla The C atoms are described by a Troullier Martinsta 
pseudo-potential with s non-locality. The wave-functions are expanded in a plane wave basis set with a kinetic 
energy cut-off of 40 Ry. The dynamical matrix is obtained by displacing each atom and computing the resulting 
forces. The Brillouin zone is sampled with the T k-point only. We compute the intensities of the resonant Raman 
spectra using Eqs. (|l2|), (|l4|), and (|l5|). The evaluation of the intensity at a resonant frequency requires a very fine 
Brillouin zone k-point mesh, for which a fully ab-initio calculation is not affordable. Therefore we evaluate the function 
k Ul .u(R, R', ^L)r-^nd the electronic susceptibility Xu L ,u(R, <^l) within a TB approximation. We use a TB model with 
s and p orbitals £3 We verified that this TB Hamiltonian well reproduces the ab-initio random-phase-approximation 
(RPA) dielectric constant of diamondEa for laser frequencies smaller than 10 eV, which is the range of energies used in 
this work for the calculation of the Raman spectra. The correct description of the dielectric constant for frequencies 
larger than 10 eV would require the inclusion in the TB model of d or s* orbitals. The Brillouin zone of the 64 atoms 
supercell is sampled with 108 special k-points. We use a value of 0.05 eV for the constant 7 in Eq. ([ti|). The derivative 
in Eq. ( |14[) is computed by numerical differentiation. 

In Fig. |2j we present our theoretical vibrational density of state (DOS). The total DOS is decomposed in its partial 
contributions coming from the sp 2 and sp 3 hybridized carbons and from the stretching and bending modes. Ej The 
sp 2 vibrations cover the entire range of frequencies of the total DOS, whereas the sp 3 vibrations form a large peak 
centered at 1000 cm -1 , which dies over 1400 cm -1 . Thus, the modes above 1400 cm -1 involve only sp 2 carbons. The 
decomposition between bending and stretching contributions shows two separated peaks, centered at 600 cm -1 and 
1000 cm -1 , respectively. Above 1000 cm' 1 the stretching modes contribute alone to the DOS. Similar theoretical 
results have been presented in Ref. for a model of hydrogenated amorphous carbon. 

Before discussing the results for the Raman spectra, we present in Fig. |^ the real and imaginary part of the 
electronic dielectric constant, e(cji) = 1 + 4ttx(^l)- Our theoretical value of 5.8 for the static dielectric constant is in 
good agreement with the value of 6.2, measured experimentally for a tetrahedral amorphous carbon sample.Eil Notice 
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that, for u> > 5 eV, the function c(u>l) is quite smooth, whereas, under 5 eV, three resonant peaks at 1.3, 3.0 and 3.9 
eV show up. The presence of such peaks, which involve transitions between localized tt states, is an artifact due to 
the limited size of our model which contains only 9 7r-bonds. This kind of features should disappear in larger models. 
The electronic structure has a gap of about 0.9 eV, above this value the incident light can be absorbed by the sample, 
as shown by the non zero value of the imaginary part of the dielectric constant. Thus, Raman spectra both in the 
visible and ultra-violet (UV) regions are collected under resonant conditions. 

Fig. H and || present the calculated Raman spectra of amorphous carbon for an incident laser beam of 1.8 eV and 
4.3 eV, respectively. The laser beam energies are chosen to be away from the three resonant peaks of e(wi). For each 
spectrum, we show a decomposition-between sp 2 , sp , bending, and stretching vibrations. We also plot the overlap 
between the different contributionsEj. The decompositions are meaningful only if the overlap is small. Again, the 
details of the shape of the Raman spectra coming from sp 2 atoms are affected by the limited number of sp 2 -sp 2 bonds 
in our model. For example the division of the contribution above 1300 cm -1 in two peaks is an artifact due to the lack 
of statistic. In good agreement with the experimental observationp the visible Raman spectrum is essentially due to 
sp 2 carbons which give a large contribution above 1200 cm" 1 , centered at 1500 cm -1 . Only a little sp 3 contribution, 
centered at 1000 cm" 1 , can be noticed. In the UV Raman spectrum, on the contrary, a contribution of sp 3 hybridized 
atoms can clearly be noticed through the large peak centered at 1000 cm" 1 . Due to a larger number of sp 3 hybridized 
atoms Wrfur model, the shape of this peak is more reliable, and its position and shape closely match the experimental 
results.Bel A peak centered at 1600 cm -1 , due to sp 2 sites, can still be observed, as in experiments. Both visible 
and UV Raman spectra are in a large majority due to the stretching vibrations, only a little contribution of bending 
vibrations can be noticed at about 600 cm" 1 , which has sometimes been observed in the experimental spectra (see 
e.g. Fig. 1 of Ref. |). 

For a given laser pulsation, ojl, we compute the integrated Raman Stokes intensities as the total area of the Stokes 
spectra. The total integrated Raman intensity and its decomposition in sp 2 , sp 3 , bending, and stretching contributions 
are presented in Figs. [| and ^ as a function of wl- These data show that whatever the energy of the incident light, 
the intensity of the Raman spectrum is always dominated by stretching modes, whose contribution is always at least 
an order of magnitude larger than that of bending modes. Regarding the relative decomposition in the sp 2 and sp 3 
contributions, Fig. || and ^ clearly show an inversion at about 5 eV. Under 4 eV almost 90% of the intensity is due to 
sp 2 carbons, whereas over 6 eV almost 90% of the intensity is due to sp 3 carbons. Notice that this inversion happens 
very quickly, thus, in an experimental measurement a little increase in the energy of the incident light in the UV 
region can lead to a major modification of the Raman spectra. 

Finally, in Fig. ||, we compare the total Raman intensities computed using our method and the Placzek's approxima- 
tion. The Raman spectra computed with the Placzek's approximation show stronger singularities in correspondence 
to the sp 2 — sp 2 electronic transitions. Instead, the intensity of the spectra computed with our approach is similar for 
laser frequency in resonance with sp 2 — sp 2 transitions and with sp 3 — sp 3 transitions. We did not find any experi- 
mental data for the total intensity as a function of the laser frequency, but according to the discussion of Section || 
the use of our approximation is more justified in the resonant case. 

IV. CONCLUSION 

In this paper, we have presented a practical method to compute the vibrational resonant Raman spectra in solids. 
Our approach is justified if the electronic excitations are delocalized. We use our method to compute the vibrational 
resonant Raman spectra of tetrahedral amorphous carbon. The computed spectra are in good agreement with the 
experimental ones measured with visible and UV lasers. In particular, the theoretical visible Raman spectra present a 
broad feature between 1200 and 1700 cm" 1 associated with sp 2 stretching as found in experiments. In the theoretical 
UV Raman spectra a second broad peak around 1000 cm" 1 associated to sp 3 stretching shows up. The location and 
the shape of this peak agree very well with the experimental UV spectra. We have analyzed in detail the evolution 
of the Raman spectra as a function of the laser frequency. We have shown that, at any frequency, the spectra are 
dominated by the stretching vibrations. We have identified a very rapid inversion in the relative Raman intensities 
of the sp 2 and sp 3 sites with the frequency of the incident laser beam. In particular, the spectra are dominated by 
sp 2 atoms below 4 eV and by sp 3 atoms above 6 eV. According to our results, it would be interesting to collect a UV 
spectrum of tetrahedral amorphous carbon with a laser frequency larger than those used in the actual experiments, 
to further emphasize the contribution of the tetrahedral sp 3 carbon sites. 

We acknowledge D. A. Drabold for providing us with the model of tetrahedral amorphous carbon. We thank Dr. 
M. Marangolo for a critical reading of the manuscript. The calculations have been performed at the Idris computer 
center of the CNRS. 



5 



1 A. C. Ferrari and J. Robertson, Phys. Rev. B 61, 14095 (2000), and references therein. 

2 Q. Wang, D. D. Allred, and J. Gonzalez-Hernandez, Phys. Rev. B 47, 6119 (1993). 

3 V. I. Merkulov et al, Phys. Rev. Letter 78, 4869 (1997). 
4 K. W. R. Gilkes et al, Appl. Phys. Lett. 70, 1980 (1997). 

5 P. Briiesch, Phonons, theory and experiments II : experiments and interpretation of experimental results (Springer- Verlag, 
Berlin Heidelberg New York, 1986). 

6 Here and in the following we will use atomic units, with which h = 1. 

7 G. Placzek, in Handbuch der Radiologie, edited by E. Marx (Akademische Verlagsgesellschaft, Leipzig, 1934), Vol. 6, p. 209. 

8 D. Porezag and M. R. Pederson, Phys. Rev. B 54, 7830 (1996). 

9 P. Giannozzi and S. Baroni, J. Chem. Phys. 100, 8537 (1994). 
10 P. V. Santos et al., Phys. Rev. B 52, 12158 (1995). 

1 R. Loudon, The Quantum Theory of Light (Clarendon Press ; Oxford University Press, New York, 1983). 

12 D. A. Drabold, P. A. Fedders, and P. Stumm, Phys. Rev. B 49, 16415 (1994). 

13 F. Mauri, B. G. Pfrommer, and S. G. Louie, Phys. Rev. Lett. 79, 2340 (1997). 

14 N. Marks et al, Phys. Rev. B 54, 9703 (1996). 

15 J. Hutter et al, CPMD Version 3.3.5, MPI fur Festkorperforschung and IBM Research Laboratory, 1990-1998. 

16 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 (1991). 

17 C. H. Xu, C. Z. Wang, C. T. Chan, and K. M. Ho, J. Phys. Condens. Matter 4, 6047 (1992). 

18 L. X. Benedict, E. L. Shirley, and R. B. Bohn, Phys. Rev. B 57, R9385 (1998). 

19 To perform the projection we have defined for each bond a 'stretching' vector in the space of the 3-/V a tom displacements. 
The components of each vector involve the displacement of two atoms in the direction of the bond and with opposite 
orientations. We use these vectors as a (non-orthonormal) basis of the stretching subspace. We define the bending subspace 
as the complement of the stretching subspace. 

20 F. Mauri and A. D. Corso, Appl. Phys. Lett. 75, 644 (1999). 

21 Z. Y. Chen and J. P. Zhao, J. Appl. Phys. 87, 4268 (2000). 

22 In the expression of intensityJEq. (§, A = fetching + ^bendmg where the contributions fetching and ^bending are 
obtained substituting in Eq. ( |14| ) the vibrational eigenvector q s with its projections in the stretching and bending subspaces, 
respectively. The intensity is then proportional to |A| 2 = |^f rctchin s| 2 + |^endtog|2 + 2 ^stretchi^benxli n g_ The first twQ 
terms correspond to the contributions of stretching and bending modes, respectively, whereas the last term corresponds to 
the overlap term. A similar repartition is used to define the sp 3 and sp 2 contributions. 




FIG. 1. A ball and stick model of the sp 2 carbon atoms in the unit cell. A bond is drawn when the distance between two 
atoms is less than 1.9 A. The sp 2 atoms are arranged in chains, no aromatic ring are present in the model. 
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FIG. 2. The total vibrational density of states of tetrahedral amorphous carbon is decomposed in the contributions coming 
from sp 2 and sp 3 carbon sites (upper panel) and in the contributions coming from bending and stretching modes (lower panel) . 
To obtain continuous curves in this figure and in Figs. ^| and ^ we use a Gaussian smearing with a = 50 cm -1 . 
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FIG. 3. Real and imaginary parts of the electronic dielectric constant of tetrahedral amorphous carbon, e(ojl), as a function 
of lo-l. 
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FIG. 4. Visible Raman spectrum of tetrahedral amorphous carbon computed for ul = 1.8 eV. The Raman intensity is 
plotted as a function of the Stokes shift, (wl — u)). The total Raman spectrum is decomposed in the contributions coming from 
sp 2 , sp 3 atoms, and their overlap (upper panel) and in the contributions coming from stretching and bending vibrations and 
their overlap (lower panel). Note that the presence of two peaks for energies larger than 1200 cm -1 is an artifact due to the 
small total number of sp 2 sites in our theoretical model. 
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FIG. 5. UV Raman spectrum of tetrahedral amorphous carbon computed for wl = 4.3 eV. The Raman intensity is plotted 
as a function of the Stokes shift, (oil — oj). The total Raman spectrum is decomposed in the contributions coming from sp 2 , 
sp 3 atoms, and their overlap (upper panel) and in the contributions coming from stretching and bending vibrations and their 
overlap (lower panel). 
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FIG. 6. Absolute integrated intensities of the Raman spectra of tetrahedral amorphous carbon plotted as a function of the 
energy of incident laser beam, Ul. In this figure, we present the absolute intensity divided by to 4 ' (see Eq. (p^)) of the 64 atoms 
model using a logarithmic scale for the intensity axis. The total integrated intensity is decomposed in the contributions coming 
from sp 2 and sp 3 atoms (upper panel) and in the contributions coming from stretching and bending vibrations (lower panel) . 
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FIG. 7. Fractional integrated intensities of the Raman spectra, plotted as a function of the energy of incident laser beam, 
u^l. Here, we present the fraction of the sp 2 and sp 3 intensity per site. For each value of wl, we compute the intensities 
per site, i sp3 (wL) and i sp2 {u-L), by dividing the absolute sp 3 and sp 2 intensities presented in Fig. ^[ by the number of sp 3 
and sp 2 atoms in the model, respectively. In this figure we plot the fractional intensities, i Bp3 (wl) /[i sp3 (ojl) + i sp2 (u>-L)] and 

^ 2 K)/[i sp3 K) + » sp2 K)]. 
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FIG. 8. The total integrated intensity of Raman spectra divided by uj a (see Eq. (|l2|)) computed with Placzek's approximation 
and with the approximation developed in this paper, as a function of the incident laser beam energy. We present the intensity 
in a logarithmic scale. 
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